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ABSTRACT 

The merger of a binary system composed of a black hole and a neutron star may leave behind a torus of hot, 
dense matter orbiting around the black hole. While numerical-relativity simulations are necessary to simulate 
this process accurately, they are also computationally expensive and unable at present to cover the large space 
of possible parameters, which include the relative mass ratio, the stellar compactness, and the black hole spin. 
To mitigate this and provide a first reasonable coverage of the space of parameters, we have developed a 
method for estimating the mass of the remnant torus from black hole-neutron star mergers. The toy model 
makes use of an improved relativistic affine model to describe the tidal deformations of an extended tri-axial 
ellipsoid orbiting around a Kerr black hole and measures the mass of the remnant torus by considering which 
of the fluid particles composing the star are on bound orbits at the time of the tidal disruption. We tune the 
toy model by using the results of fully general-relativistic simulations obtaining relative precisions of a few 
percent and use it to extensively investigate the space of parameters. In this way we find that the torus mass 
is largest for systems with highly spinning black holes, small stellar compactnesses, and large mass ratios. As 
an example, tori as massive as Af & tm — 1-33 AIq can be produced for a very extended star with compactness 
C ~ 0.1 inspiralling around a black hole with dimensionless spin a = 0.85 and mass ratio q ~ 0.3. However, 
for a more astrophysically reasonable mass ratio q ~ 0.14 and a canonical value of the stellar compactness 
C ~ 0.145, the toy model sets a considerably smaller upper limit of Afbtor ^ 0.34 Mq. 

Subject headings: gamma rays: bursts - black hole physics - stars: neutron 



L INTRODUCTION 

The most widely accepted scenario to explain the phe- 
nomenology associated with short-hard gamma-ray bursts 
(SGRBs) involves the merger of either black h ole (BH) neu - 
tron star (NS) systems or of binary NS systems (Nakar 2007| l. 
In either case, the remnant consists of a BH with negligible 
baryon contamination along its polar symmetry axis and of 
a hot, massive accretion torus surrounding it, which releases 
energy as it accretes onto the black hole, typically in the form 
of a relativistic jet. With these fundamental ingredients of the 
standard SGRB model, an intense neutrino flux is emitted as 
the torus accretes onto the BH, triggering a high-entropy gas 
outflow off the surface of the accretion torus, i.e. "neutrino 
wind". At the same time, energy deposition by vD annihi- 
. . lation in the baryon-free funnel around the BH rotation axis 
>- powers relativistically expanding e* jets, which can give rise 



tions have made considerable progress in t h e last few years 
(see, f o r instance [OechsUn & Janka| (|2007|l; [Anderson et al. 



(|2008l); 'Baiotti e t al.|(|2008 ');'Yamamoto et al.|(|200 8); "Ander- 



son et a l. (2008);iLiuetal. (20 08); Giaco mazzo et al. (2009ji^ 
Rezzolla et al. (2010); Bauswe in et al.|(|2 010) for recent stud- 
ies of NS-N S binaries or Kluzmak & Lee (1999); Rosswog 
(|2005|);|Loffl er et al. (2006); Etienne et al. (2008); Shibata S. 
"Taniguchi'C2008); Duez et al. (2010); Etienne et al. (20091; 
Duez et al. (2008); Shibata et al. (2009); Chawla et al. (201()| 
for corresponding work on BH-NS binaries). Despite the fact 
that this type of simulations is now possible, they remain nev- 
ertheless both challenging and computationally intensive. Nu- 
merical simulations of NS-NS mergers have now reached a 
rather high level of accuracy (see Baiotti et al. (2009) and 
the discussion in the Appendix of JRgzzolla et al.| ( |2010p ), and 
different codes have been shown to yield results th at agree 



to the observed GRB emission. Other burst mechanisms have, 
of course, been proposed; since these principally involve mag- 
netically launched jets and since we do not address magnetic 
fields in this paper, we have briefly summarised only the burst 
mechanism powere d by neutrino annihilat i on; th e interested 
reader may refer to Lee & Ramirez-Ruiz ( 2007| l for a thor- 
ough review. 

The simulation of these events "ab initio" requires an 
adequate description of general relativity, relativistic (mag- 
neto)hydrodynamics, and a proper microphysical equation of 
state. Typically, the only way to model these systems accu- 
rately is to resort to numerical-relativity simulations, solving 
consistently both the Einstein equations and those of relativis- 
tic hydrodynamics or magnetohydrodynamics. These simula- 
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to 10% (at worse) when using the same initial data ( [Baiotti 
[et al.|[20T0| l. However, the situation is much more uncertain 
in the case of BH-NS binaries, for which no direct compari- 
son among different codes has been made yet and the results 
of the simulations from different codes are sometimes not in 
agreement. As an example, the merger of the same binary 
with mass ratio 1/3 yi elds a torus w i th bar yon mass which 



is ~ 4% of the NS in Etienne et al. (2009) and < 0.001% 



in [Shibata et al.[ ( [2009| l. As a result, no reliable knowledge 
is available at the moment on how the mass of the torus de- 
pends on the most important parameters of the system: the 
mass ratio, the stellar compactness, and the BH spin. 

These problems, along with the need of a better under- 
standing of the tidal disruption process, have pushed the par- 
allel development of ps eudo-Newtonian BH-NS calculations 
- e.g. 'Ruffert & Janka ( 2010J ) use the Paczyriski-Wiita phe- 
nomenological potential to mimic the innermost stable circu- 
lar orbit (IS CO) of the BH in a Newtonian setting - and of 



se mi-analyt i cal ap proaches to the problem. Regarding the lat- 
ter, Shibata ( 1996| l, for instance, described the necessary con- 
ditions for the production of an accretion torus of appreciable 
size by requiring that the NS disruption occurs at a tidal radius 
rtidc that is larger than the ISCO of the BH r-jsco- Unfortu- 
nately, Shibata (199 6) did not predict the mass of the result- 
ing accretion torus except to assume that it vanishes when the 
radius of tidal disruption is less than that of the ISCO. A par- 
allel systematic study has been pursued recently to exploit the 
relativistic "affine-m odel" and desc ribe the prop erties of th e 
tidally deformed NS ( jFerrari et alTT2009 ; Ferrari et al |2010| l. 
Also in this case, however, the study has concentrated on the 
evolution of stationary configuration and has not made any 
prediction on the final outcome of the merger. 

In this paper we attempt to bridge the gap between intensive 
numerical simulations and semi-analytical studies by estab- 
lishing a way to estimate the mass of the resulting torus. We 
do this by taking the concept of the tidal disruption to its log- 
ical extreme. In other words, we model the NS in the binary 
as a relativistic tri-axial ellipsoid which is tidally distorted as 
it orbits in the tidal field of a rotating BH. When the tidal- 
disruption radius is reached, however, we assume the star to 
be composed of a system of non-interacting "fluid particles" 
which move on the corresponding geodesies. We therefore 
compute the mass of the torus as the integral of the masses of 
the particles which do not fall into the BH. This clearly rep- 
resents only a "toy model" for the complex dynamics of the 
merger process, but we show that, with a suitable tuning, it al- 
lows us to reproduce with good precision the large majority of 
the results obtained so far from more accurate but also consid- 
erably more expensive numerical relativity calculations. Most 
importantly, however, it provides a simple tool to better un- 
derstand the complex dynamics of the tidal disruption and to 
cover at once the full space of parameters. 

The structure of the paper is the following one. In Sect. |2] 
we describe the particular tidal model we use and then how we 
estimate the mass of the accretion torus. In Sect. |3] we show 
that by tuning the free parameter in our model we can repro- 
duce results obtained within fully general-relativistic simula- 
tions, thus proving that the tool we build is solid. In Sect.|4]we 
present the results of our estimates, leaving an intuitive inter- 
pretation of the results and the conclusive overview to Sects. [5] 
and l6] respectively. 

2. METHOD 

To model the behaviour of the NS during the final stages of 
the inspiral of the mixed binary and before it merges with the 
BH, we use an improved version of t he affine model which 
is presented in detail in Ferrari et al.| (12009 ). An important 
difference with respect to that work is that we do not consider 
the prescriptions of the quasi-equilibrium approximation and, 
rather, follow the dynamics of the NS until it is disrupted by 
the BH tidal field. Furthermore, as mentioned in the Intro- 
duction, in addition to treating the NS as a tri-axial ellipsoid, 
we also decompose it into a large number of representative 
"fluid particles", whose kinematic properties will be used to 
study the motion of the NS matter after the tidal disruption. 
Our toy model is therefore composed of three logical parts: 
(i) the evolution of the NS deformation as it inspirals towards 
the BH; (ii) the modelling of the tidal disruption; (Hi) the cal- 
culation of the mass building-up the torus. Each of these parts 
will be discussed separately in the remainder of this Section. 

2.1. Neutron star deformation - The affine model 



The idea of modelling stars as ellipsoids has a long his- 
tory and a thorough analysis of incompressible ellipsoidal 
figures of equilibrium was performed by Chandrasekhar in 
1969 (Chandrasekhar 1969). When modelHng the NS defor- 
mation, we are addressing what is known as the compressible 
Roche-Riemann problem, in which one studies the behaviour 
of a compressible ellipsoid with uniform vorticity parallel to 
its rotation axis, orbiting a point mass or a rigid sphere. More 
specifically, we will be using an improved version of the affine 
model, which was developed in the '80s by Carter, Luminet 
and Marck to describe the encounters between a BH and a 
Newtonian s tar (Carter & Luminet 1982, 1983 1985;Luminet] 
& Marckll"985, ,Luminet & Carter 1986 ) and then applied to 
BH-N S bmaries at the very end of the '90s (Wi ggins & Lai| 
2000 1. More recently, the Newtonian treatment of the star was 
upgraded to achi eve a better description of the NS in m ixed 
compact binaries ( [Ferrari e t al. 2009; Fer rari et al.|2 010). The 
essential features and assumptions of the improved affine ap- 
proach used here can be summarized as follows: 



• the equilibrium structure of the NS is determined by the 
TOV equations, while its dynamical behaviour is gov- 
erned by Newtonian hydrodynamics improve d by the 
use of an effec tive relativistic scalar potential ( |Rampp| 
& Janka|2002l l; 



• the NS centre of mass moves in the tidal field of a 
Kerr BH along a simple inspiralling equatorial orbit 
{cf. eq. (21 1] and each point of the orbit is associated 
with a BHtimelike circular geodesic having the same 
radius; 



• throughout the inspiral, the NS remains a Riemann S- 
type ellipsoid, i.e. its spin and vorticity are parallel and 
their ratio is constant (see ,Chandrasekhar (1969^); 



• tidal effects on the orbital motion and the perturbation 
that the star induces on the BH are neglected. 



For completeness, we next review the mathematical formu- 
lation of the affine model used here by writing the equations 
governing the NS deformations in the principal frame, i.e. the 
frame associated with the principal axes of the stellar ellip- 
soid. In this frame, the fluid variables of the affine model are 
five: the three principal axes of the stellar ellipsoid oi, 02, and 
as, the angular frequency of the internal fluid motion A and 
the star spin Q. measured in the parallel-transported frame as- 
sociated with the centre of mass of the star ( Marck 1983 ). The 
axis 03 is perpendicular to the orbital plane, while Oi and 02 
are perpendicular to a-^ (see Fig. [Til. In the Newtonian limit 
(AfsH ^ r, where Mbh is the BH mass and r is the Boyer- 
Lindquist radial coordinate), ai and 02 belong to the orbital 
plane. The dynamics of four of the five fluid variables is then 



governed by the following set of equations 
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where the dots denote derivatives with respect to the proper 
time at the stellar centre r, the Cij 's denote the components 
of the BH tidal tensor in the principal frame, i?Ns is the NS 
radius, and the index symbols Ai are defined as^ 
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The effective relativistic self-gravity potential for the isolated 

NS in spherical equilibrium equilibrium V and the scalar 
quadrupole moment for the isolated NS in spherical equilib- 
rium equilibrium M are given by"* 

r-RNs [e(f) +p(f)][mTov(^) +47rr''p(f)] ^3^ 
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is the spin angular momentum of the star. The fifth fluid vari- 
able may be expressed in terms of 
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which is proportional to the circulation in the locally nonro- 
tating inertial frame. We note that because we work in ab- 
sence of viscosity, the circulation of the fluid is conserved, 
i.e. C = 0. For simplicity, but also because this is the assump- 
tion made by all numerical-relativity simulations to date, we 



The A,'s are related to the diraensionless index coefficients defined 



Chandrasekhar 



^„llalnJlao^l^Jllll jl969) by the simple dimensional rescaling Ai = R^^Ai. 

essence, they express the derivative of the self-gravity of the deformed star 
with respect to its i-th axis. 

* Hats ( " ) denote quantities calculated for an isolated nonrotating NS at 
equilibiium. 
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Fig. 1. — Cartoon of the toy model. Indicated are the tidal radius r^^^^, 
the ISCO risco: two of the principal axes ai,a2, the principal frame and the 
parallel-transported frame. Note that for simplicity we set = ^. 



set C = initially (the NS fluid is thus irrotational) so that 
( [TO] i reduces to 
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The components of the tidal tensor for a Kerr spacetime, 
expressed in the NS principal frame, are 
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where the angle (p (which is related to il hy <j) = il) is the an- 
gle that brings the parallel-transported fr ame into the p rincipal 
frame by a rotation around the 03 axis ( |Marck|1983| l. Simi- 
larly, ^ is an angle that governs the rotation of the parallel- 
transported tetrad frame in order to preserve the parallel trans- 
port of its basis vectors. Stated differently, the difference be- 
tween ^ and (p represents the lag angle between the principal 
frame and the parallel transported one, i.e. (j>ii^g = 'i' — (f>, and 
thus measures how much the star is "lagging behind" in its 
orbit around the BH. If (f)\^g — and thus 'i' — (f), then the 
largest semi-major axis of the star is always pointing towards 
the BH. A schematic diagram of the BH-NS binary and of the 
relevant quantities discussed so far is shown in Fig. [T] 

The constant K appearing in the tidal-tensor compo- 
nents ([T2l)-([T5|) is a combination of the energy E and the z- 
orbital angular momentum per unit mass of the star Lz 
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where a = 
geodesies. 



J/M is the spin of the BH and, for circular 



E: 



2rMn 



ay/rMs 



2ay/rMBH + a^) 






v/r-2 - 3r-MBH + 2ay/iMi 



(17) 
(18) 



BH 



Note that because the tidal-tensor components cn and C22 
have a different sign, the forces acting on the corresponding 
semi-major axes oi and 02 also have opposite signs and thus 
lead to a stretching of oi and a to compression of 02- 

For simplicity, and to obtain a better agreement with the 
results of the numerical-relativity simulations, we will con- 
sider hereafter (/)iag — and thus = ^. Furthermore, since 
we will consider a sequence of circular equatorial geodesies 
whose radius reduces due to the emission of gravitational ra- 
diation, the evol ution o f the angle 'J for each circular orbit is 
given by (Marck|1983|l 



and thus also 
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In order to evolve the equations of the affine model ([T]i-(|3]l 
we must select an equation of state (EOS) for the NS matter 
and specify the initial conditions and the evolution of the or- 
bit. As far as the first is concerned, the model is sufficiently 
general that any EOS co uld be used and ind eed several dif- 
ferent ones were used inj Ferrari et al.| ( |2010| l. However, be- 
cause here we want to compare with the results of numerical- 
relativity simulations and these have been performed mostly 
with a r law equation of state p — (T ~ l)pe which for the 
adiabatic process considered in this paper is equivalent to a 
poly tropic EOS p = Kp^ with F = 2 or 2.75, we will con- 
sider here just polytropes with these polytropic exponents and 
present the result of more realistic EOSs in a subsequent work. 
As for the initial conditions, we consider an initial separation 
for the binary, vq = r{t = 0) and set ^o = (/)(t = 0) = 0, 
while n, and A are automatically given by ( |20] i and ( flT) , re- 
spectively. For the NS axes, instead, we set the time deriva- 
tives on the left hand sides of ([T)-<l3j to zero and solve the sys- 
tem for oi, 02, 03 with a Newton-Raphson scheme. Of course, 
it is necessary to ensure that rp is large enough, i.e. that ini- 
tially fli ~ a2 — as ~ i?NS and Js ~ (as well as C = 0), 
and that the final results are unchanged if one chooses a larger 
tq. In other words, att = Q the NS must almost be at spher- 
ical equilibrium and the calculations must therefore be inde- 
pendent of the specific choice made for vq. 

For the time evolution of the orbital separation r, we con- 



sider a very simp le circular equatorial adiabatic inspiral ( Mis- 
|ner et alI1|1973| l, which accounts therefore for the radiative 
losses of two point masses at a 2.5 post-Newtonian (PN) ap- 
proximation 



r{t) = ro 1 
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where 
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is the inspiral time and A/ns is the gravitational mass of the 
NS. When generating the orbit, we evolve the orbit angle 9 
according to the Kerr spacetime equation 



de _ 1 

dt a + ^r'VA/B 



and we make the approximation of setting 
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It is important to remark that our goal is that of computing the 
mass of the torus produced by the tidal disruption and not that 
of providing an accurate description of the binary inspiral. In 
this sense, using a lower-order PN description of the orbit is 
very reasonable as the dynamics we are most interested in take 
place when the presently available PN models are no longer 
accurate. 

Once the orbit is determined, we may integrate the affine- 
model equations, terminating the evolution when the ratio of 
the semi-major axes reaches a critical value (a2/ai)ci-it. This 
quantity cannot be determined a priori and is effectively a free 
parameter in our toy model. However, it may be tuned by 
comparing the results of the toy model with those of the nu- 
merical simulations and the way we do this will be discussed 
in the next Section. We thus define the tidal disruption radius 
r,i(ie as the orbital separation at which (a2/ai) — (a2/ai)crit. 

A final quantity which is relevant to introduce and that may 
be useful to interpret the results of the toy m odel is the ISCO, 
which, for a generic Kerr BH is given by (Bar deen et al.|1972) l 

nsco = A/bh{3 + Z2T [(3 - ^i)(3 + Zi+ 2Z2)]i/2| ^ 
Zi = 1 + (1 - aV^B^H)'/' , 

X [(1 + a/AfBH)'/' + (1 - S/A//bh)'/'] , 
Z2 = {'SdyMin + Zfy/\ (25) 

where the upper/lower sign holds for co-rotating/counter- 
rotating orbits. In general, the ISCO is inside the tidal radius, 
i.e. risco < ?'tide, but there are situations in which the oppo- 
site is true and this is the case, for instance, when considering 
binary systems with very small mass ratios or stars with very 
large compactness. In these cases too, we follow the evolution 
of the axis ratio and "disrupt" the NS inside the ISCO as soon 
as the critical value is reached. 

2.2. Neutron star disruption 

As mentioned above, when the affine-model evolution of 
the mixed binary leads to the tidal disruption of the NS, we 
fragment the NS into fiducial fluid elements that would be 
representative of the motion of the NS matter. The first step in 
our strategy consists therefore in switching from the five fluid 
variables of the affine-model formulation to a description of 
the (disrupted) NS fluid as a set of test particles, each one of 
which possesses a mass and a 4-velocity. In practice, at dis- 
ruption we build a fine grid adapted to the ellipsoidal shape 
of the star and divide the star into a collection of fluid ele- 
ments. In the principal frame, the centre of each fluid element 
is identified by a 3-vector x and we calculate the mass of the 
corresponding fluid cell by multiplying the mass density at its 
centre by the volume of the fluid element. Moreover, we may 
associate to the centre of mass of each cell a 3-velocity, which, 
in the principal-axes frame, is given by (Chandrasekhar^ 1 969) 
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is the spin velocity (/.e. the speed of the fluid due to its rota- 
tion), and 



Ue = — xiei 
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is the eflipsoid expansion/contraction velocity, e^ being the 
unit vectors along the ellipsoid principal axis a.i. The coor- 
dinate Xi along the i-th axis runs from — a^ to a^. With a 
rotation of an angle (f) around 03, we switch from position 
3-vectors in the principal frame to position 3-vectors in the 
parallel-transported tetrad associated with the NS centre of 
mass, where the 3-velocity u becomes u + Cl x x. In this 
parallel-transported tetrad, we determine the time component 
of each 4- velocity vector by exploiting the normalization con- 
dition u^"^U(^a) = ^1 and by recalling that the time compo- 
nent of the position vectors is simply x^"-* = 0. Finally, we 
express all the 4-position and the 4-velocity vectors in Boyer- 
Lindqui st coordinate s by applying the transformation laws de- 
rived in'Marck'(1983) and summarised in Appendix [A] 

The procedure described above provides a complete de- 
scription, in Boyer-Lindquist coordinates, of the kinematic 
properties of fluid parcels as point particles freely-falling in 
a Kerr spacetime; this is what is needed to then estimate the 
torus mass. 

2.3. Torus mass estimation 

When looking carefully in numerical-relativity simulations 
at the dynamics of the NS after it is disrupted, it is quite strik- 
ing to note how much the different parts of the star seem to 
behave like independent freely falling particles: the gravity 
of the black hole alone does seem to represent the dominant 
force at this stage of the evolution. In view of this observa- 
tion, when the NS is tidally disrupted and split into fiducial 
fluid elements of which we know the mass and the 4-velocity, 
we assume that the pressure gradients across neighbouring el- 
ements and the self-gravity of the system play little role, and 
hence that the fluid elements behave as independent collision- 
less fluid particles. As such, after the disruption the NS is 
approximated as an ensemble of about 3.1 x 10"* fluid parti- 
cles which have a complex distribution of energy and angular 
momenta, but are in free-fall towards the BH. 

Using the 4-velocity, of each particle we compute the cor- 
responding conserved quantities by inverting the relations 
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and e, £z, and Q represent the energy, angular momentum and 
Carter's constant of motion, respectively, all normalized to the 
mass of the particle. Note that it is necessary to use these 
general equations instead of equations ^Tl\ and ( [T8] l as the 
majority of the particles no longer follows circular equatorial 
geodesies. 

As mentioned previously, we identify the mass of the rem- 
nant torus with the sum of the masses of the bound particles 
and we make use of ([30| to determine whether a given parti- 
cle is bound or not. Noting that a turning point occurs when 
[dr/drY passes through zero and since the only influence of 
6 is to decrease the overall magnitude of {dr/drY but not 
to change its sign, we only consider, without loss of general- 
ity, the case 6 — 7r/2. We then use root-finding techniques 
for each particle and consider bound those particles for which 
[dr/drY < at a radial position rxp outside the event hori- 
zon Teh, such that r^H < ?'tp < ?'tide (note that [dr/dry is 
always greater than at r = rtide) and simultaneously satisfy 
e < 1. This final condition merely states that the gravitational 
binding energy has compensated the kinetic energy such that 
the total energy of the particle is less than the rest mass of the 
particle at infinity. 

Once the NS is tidally disrupted, the calculation of the torus 
mass, which is initially set to be Afb.tor — ^j^b.NS' is done as 
follows. 

1. For each fluid particle we verify whether it is bound 
or not. In this latter case we assume the particle will 
accrete onto the BH^ . 

2. The composite mass of the accreted particles is added 
to the mass of the BH and the mass of the torus is de- 
creased by the corresponding amount. 

3. We reconsider the remaining particles and verify if they 
are still bound or if they would now accrete onto the 
new and more massive BH. 

This procedure is repeated until there are no more particles 
that would accrete onto the BH or, equivalently, until the rel- 
ative change in the mass of the torus is less than one part in 
one million. 

In addition to a change in the mass of the BH we have also 
experimented with changing the spin of the black hole as a 
result of the angular momentum accreted with the particles. 
However, the results in this case are much less robust (the 
mass of the torus is not a monotonic function of the parame- 
ters) and this is probably due to the more complex dependence 
of the geodesic motion on the spin of the BH, which conflicts 
with the approximations made here. As a result, we keep the 
BH spin to be the same as the initial one and it is reassur- 
ing that this does not spoil the very good agreement with the 
numerical simulations. 

3. TUNING AND VALIDATION OF THE MODEL 

In the affine-model approach base d on a quasi-equ i librium 
approximation and discussed in [Ferrari et al.| ( |2009[ ), 
the disruption radius is identifieH By tEe condition 
[d{a2/ai)/dr]~^ = 0, i.e. as the radial separation at 
which the axis ratio diverges. Although this singular limit is 

^ We note that all numerical simulations suggest that the amount of matter 
leaving the central gravitational potential, i. e. that are unbounded but do not 



fall onto the BH, is extremely small and can thus be neglected here (see jRez-] 
|zollaetal.|pOTO) ). 



TABLE 1 

Comparison between the remnant torus mass predictions of fully general-relativistic simulations and of our model with the 

critical value of a^jai tuned to 0.44. in the four sections of the table, we examine the results recently provided in (from top to 

bottom) TONITA ET AL. (2010), DUEZ ET AL. (2010), ETIENNE ET AL. (2009), AND SHIBATA FT AL. (2009). THE FIRST FOUR COLUMNS OF THE 
TABLE ARE THE PARAMETERS OF EACH BH-NS BINARY, i.e. THE ADIABATIC INDEX V OF THE NS EOS, THE NS COMPACTNESS C, THE MASS RATIO q 
AND THE DIMENSIONLESS BH SPIN a. THE FOLLOWING THREE COLUMNS PROVIDE THE REMNANT TORUS MASSES A/tor OBTAINED FROM OUR MODEL 
(LABELLED "toy model"), THOSE OBTAINED FROM FULLY GENERAL-RELATIVISTIC CALCULATIONS (LABELLED "simulations"), BOTH GIVEN IN UNITS 

OF THE NS BARYONIC MASS Mt, nS. AND THE RELATIVE ERROR. 



Ref. 


EOS 


C 


1 


a 


.'^^b.tor/Afb.NS 


A4,tor/A^(j,NS 


error 




(F) 








{toy model) 


(simulations) 


(%) 


Tonita et al.'feOlO'i 


2.00 


0.100 


1/5 


0.00 


0.17 


0.17 





Tonitaetal. (20101 


2.00 


0.125 


1/5 


0.00 


0.06 


0.06 





Tonita etal. (2010) 


2.00 


0.145 


1/5 


0.00 


< 0.01 


<0.01 





Tonitaetal. (20101 


2.00 


0.150 


1/5 


0.00 


< 0.01 


<0.01 





Duez etal. (20101 


2.00 


0.144 


1/3 


0.50 


0.08 


0.08 





Duez etal. (20101 


2.75 


0.146 


1/3 


0.50 


0.11 


0.13 


18 


Duez etal. (20101 


2.75 


0.173 


1/3 


0.50 


0.04 


0.02 


50 


Etienne et al. 1 2009 ) 


2.00 


0.145 


1/3 


0.00 


0.02 


0.04 


100 


Etienne et al. 1 2009 ) 


2.00 


0.145 


1/3 


0.75 


0.18 


0.15 


17 


pienne et al. 1 200^ 


2.00 


0.145 


1/3 


-0.50 


< 0.01 


<0.01 





Etienne et al. 1 2009 


2.00 


0.145 


1/5 


0.00 


< 0.01 


<0.01 





Shibata et al. 1 20091 


2.00 


0.145 


1/3 


0.00 


0.02 


<0.01 


100 


Shibata et al. 1 2009 1 


2.00 


0.160 


1/3 


0.00 


< 0.01 


<0.01 





Shibata et al. 1 2009 1 


2.00 


0.178 


1/3 


0.00 


< 0.01 


<0.01 





Shibata et al. 1 2009 1 


2.00 
2.00 


0.145 
0.145 


1/4 
1/5 


0.00 
0.00 


0.01 
< 0.01 


<0.01 
<0.01 


100 



Shibata etal.|j2009^ 



clearly a shortcoming of the assumption of quasi-equilibrium, 
it is not obvious how to specify the tidal radius in a way 
which is not arbitrary to some extent. To remove at least 
in part this degree of arbitrariness, we have decided to tune 
the tidal radius by carefully analyzing the results of recent 
numerical -relativity sim ulations and in p articular those 
carried out at the AEI ( Tonita et al.||2010) l, for which we 
have more direct control over the errors. When doing so, 
we realized that the critical value of the axis ratio (a2/ai)(.,it 
is a robust measure across our simulations, but also when 
comparing with the simulations published in the literature. 
Hence, we have decided to consider the critical axis ratio 
(a2/ai)ci-it as a free parameter and to identify its value as the 
one which best reproduces the numerical data available. 

More specifically, for those initial data for which numerical 
simulations have been performed, we tuned, within the toy 
model, the choice of the free parameter (a2/ai)crit so as to 
minimize the difference, with the corresponding numerical- 
relativity result, for the torus mass. As a result of this proce- 
dure we we obtain (02/01 )crit = 0.44 which is robust across 
all of the simulations and thus define the tidal radius as the or- 
bital separation at which (02/01) attains the critical value. It 
is reasonable to expect that (o2/oi)ciit will depend on the BH 
spin and on the mass ratio. Here, however, we assume that 
such dependence is weak and thus set it to be be constant. As 
we discuss below, even with this crude approximation we can 
reproduce most of the numerical results with an error which 
is below ^ 15%. As an addition note, we stress that although 
robust {i.e. a single choice fits well all of the available data), 
the masses of the tori are rather sensitive to the choice for the 
critical axis ratio. In particular, for the same binary, a change 
of ^ 2% in (o2/oi)(;rit (i.e. a change in the last significant 
figure) may lead to a change in the last significant figure of 
the estimated torus mass, and thus up to a ~ 50% change for 
cases with a very small remnant mass. This effect disappears 
if one tunes (o2/oi)ci-it with an extra significant digit. 



Before going to the details of the comparison with the nu- 
merical simulations it is worth making two remarks. The first 
one is that after having identified in the axis ratio a consis- 
tent parameter which we constrain to the second significant 
figure, we also expect that it will be further refined as new 
and more accurate results from numerical simulations become 
available. The second one has already been made in the intro- 
duction and stresses the fact that the numerical data itself does 
not show a great degree of consistency. While there are two 
cases which have been considered by more than one group, 
most of the data available refers to configurations which are 
slightly different and hence difficult to compare. Even the 
actual procedure followed to measure the mass of the tori 
differs from group to group; while most decide to measure 
the mass at a given time after the formation of the apparent 
horizon, not all groups use the same time. It would certainly 
be more reasonable if the measure was performed only when 
the mass accretio n rate has reach ed a very small and constant 
value, as done in Rezzolla et al. (2010), rather than setting a 
time which may vary from simulation to simulation. Notwith- 
standing these difficulties, it is remarkable that even for the 
same configurations (cf. the eighth and the twelfth rows in Ta- 
ble [Til, or some which are not very different (cf. the fifth and 
ninth rows in Table [Til, the numerical results yield tori whose 
masses differ considerably. Interestingly, the predictions of 
the toy model are equally distant from the numerical results 
reported in the eighth and the twelfth rows, thus suggesting 
that both simulations may be equally imprecise. 

Being a toy model, its validity is constrained to within spe- 
cific ranges of the space of parameters, which we discuss be- 
low and which allow us nevertheless to cover essentially all 
of the complete space of parameters. The first constraint on 
the range of validity comes from the mass ratio, which cannot 
be too large since the affine model assumes that the NS inspi- 
rals as a test fluid and is therefore increasingly more accurate 
the smaller the mass ratio. As a result, we will consider only 



TABLE 2 

Gravita tional-wave frequency at t he onset of tidal disruption /°^ as computed with our model (labelled "toy model") or as 

QUOTED in'Shibata & TaniguchiV2008^ (labelled "simulations"). The latter values were calculated by means of a fitting formula 

determined in Taniguchi et al. ( 2008 ) BY using quasi-equilibrium sequences of mixed binaries in circular orbits, obtained by 

solving the Einstein constraint equations in the conformal thin-sandwich decomposition and the relativistic equations of 

hydrostationary equilibrium. 



Q 


Afb.NS [Mq] 


Mns VMq} 


i?NS [km] 


(toy model) 


f2r. [kHz] 
(simulations) 


0.327 


1.400 


1.302 


13.2 


0.856 


0.855 


0.327 


1.400 


1.294 


12.0 


0.997 


0.993 


0.328 


1.400 


1.310 


14.7 


0.736 


0.738 


0.392 


1.400 


1.302 


13.2 


0.877 


0.867 


0.392 


1.400 


1.294 


12.0 


1.021 


1.010 


0.281 


1.400 


1.302 


13.2 


0.840 


0.843 



binaries with mass ratios 

0.10 < q< 0.33. 

W hile this condi t ion rem oves several of the values reported 
in IShibata et al.| (J2009), it is not at all unrealistic. We re- 
call, in fact, that the most recent estimates for the mass ac- 
creted onto the primary compact object during the common- 
envelope phase are rather low and thus the BH masses in 
close BH-NS binaries are likely to fall prim arily in values 
near AffiH — 10 M© ( Belczynski et al.||2007 l. Considering 
a canonical 1.4i\/0 NS, BH-NS systems are therefore most 
likely to come in a mass ratio that is q ~ 0.14. 
The second constraint comes from the stellar compactness, 
which cannot be too small for a relativistic compact star, nor 
too large given the test-fluid hypothesis of the affine model. 
As a result, we will consider only binaries where the NS has 

0.1 < C<0.16. 

This range covers well the one considered so far in numeri- 
cal simulations {cf. Table [Tl), but it is worth re marking that the 



( |2010| suggest a rather 
,16, which is at the edge 



recent arguments made in Ozel et al. 
high and generic compactness, C ^ 
of the range cons idered here, and t hat a standard cold EOS, 
such as the APR ( [Akmal et al.|199 8) EOS, leads on average 
to compactnesses C ~ 0.18, thus outside of the range consid- 
ered here. Future simulations in which this EOS is employed 
will help extend the range of validity in compactness of the 
toy model. 

The third constraint comes from the BH spin, which we can- 
not take as too large given that we treat the motion of the dis- 
rupted NS with geodesies and these would lead to incorrect 
results if the dimensionless spin parameter a = J /M^y^ is too 
high (e.g. the ratio A/^.tor/Affc ns ^^ 1 f or a ^- 1). As a result, 
we will consider only binaries where the BH has 

0.0 <a < 0.85. 

The fourth and final constraint comes from the mass of the 
torus, for which we need a lower limit. This is even true for 
numerical simulations, whose precision is not infinite. As a 
result, we consider the tori to have a zero baryon mass if 



Mb tor < 0.01 M, 



6,NS 



0.014 M, 



■ 



The results of the comparison are summarized in Table [T| 
where, in addition to our numerical simulations ( jTonita et al. 



2010 1, we have considered also the data reported in 



Duez et al. 



columns of the table: these are the adiabatic index V of the NS 
EOS, the NS compactness C, the binary mass ratio q, and the 
initial BH spin a, which does not change from its initial value 
in our toy model. The last three columns provide, instead, 
the baryonic torus masses A/f,_tor obtained from the toy model 
or by the fully general-relativistic calculations (both in units 
of the NS baryonic mass A/f, ns), and the relative percentage 
error. 

A rapid inspection of the Table and in particular of its last 
column clearly shows that there are four out of sixteen numer- 
ical relativity results that have rather large errors, i.e. between 
50 and 100%. Not having a clear measure of the error associ- 
ated with those simulations, it is hard to judge whether this is 
a limit of the toy model or whether this is a limit of the numer- 
ical simulations. It should be remarked, however, that these 
simulations are those which report tori masses that are close 
to the limit we consider reasonable {i.e. Af&,toi7Aff,,NS — 0.01) 
and clearly new simulations of those binaries are necessary to 
settle these differences. However, with the exception of those 
cases, the Table also reveals that the toy model can reproduce 
the remaining cases (which represent three quarters of the data 
available) with an error which is at most 18% and is virtu- 
ally for most of the cases. Considering that the numerical- 
relativity simulations in Table[T]were performed with different 
codes, different initial separations and amounts of eccentric- 
ity'', we believe that the tuning made for the toy model is both 
reasonably robust and accurate. 

To further validate the model, we use it to determine the 
frequency of the gravitational radiation emitted at the on- 
set of the tidal disruption. Results f or these f requencies 
were provided, for example, by Shibata & Taniguchi (2008) 
and Taniguchi et al. (2008), and we compare to the former m 
Table [2] and to the latt er in Fig. l2] The frequencies given in 
Taniguchi et al. ( 2008 1 are foundby using quasi-equilibrium 
sequences of mixed binaries in circular orbits, obtained by 
solving the Einstein constraint equations in the conformal 
thin-sandwich decomposition and the relativistic equations 
of hydrostationary equilibrium; a fitting formula for the fre- 
quency is also pro vided and this is used, in turn, in |Shibata & 
Taniguchi ( 2008] ). All cases considered by these authors re- 
fer to nonspmning BHs and irrotational NSs, so that, for our 
model, the gravitational-wave frequency at tidal disruption is 

* Although not often discussed, the presence of eccentricity in the initial 
dat a can lead to signifi cant changes in the mass of the torus in BH-NS merg- 
ers J^mtaeFaTpOTOl. 



PJTD| ); [E5enne et al.| ( [2009l l, and [Shibata et~aL] (|2009'). The 
parameters of each BH-NS binary are reported in the first four 
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Fig. 2. — Gravitational-wave frequency at tidal disruption shown as a func- 
tion of the (inverse of the) mass ratio q (here iJr is the polytropic length 
scale A"^' (^'"~^)). Data plotted with diamonds is produced with our toy 

• • "{2008) . 

All neutron stars have F = 2, while their baryonic mass is mdrcated by the 
colour code in the legend. This figure should be compared with the coiTe- 
sponding Fig. 1 in Ferrari et al. (2009). 



model, while data plotted with triangles is taken 



from |Taniguchi et al.| f 
jnic liiass is mdicated 



given by [cf. eq. (23 i] 



J tide 



1 Mr 



(35) 



i.e. by twice the Schwarzschild orbital frequency at the tidal 
disruption radius. We note that Fig. [2] reports the (inverse 
of the) mass ratio q as a function of TrTtr/t^*, where R^ is 
the polytropic length scale ii'i/(2r-2) Qata plotted with di- 
amonds is produced with o ur toy model, w hile data plotted 
with triangles is taken from |Taniguchi et aL| (20Q8). All neu- 
tron stars have F = 2, while their baryonic mass is indicated 
by the colour code in the legend. An inspection of Table [2] 
and of Fig. l2] shows that a very good agreement is obtained 
not only for the torus mass, but also for the gravitational- 
wave frequency. Stated differently, the assumptions that go 
into our toy model allow us to accurately capture both the or- 
bital evolution soon before the NS disruption takes place and 
the dynamics of the matter after the NS has been disru pted. 
Moreover, a com parison between our Fig.l2]and Fig. 1 in Fer 



rari et al. ( 2009|l shows that the present implementation of 



the affine model is significantly improved with respect to its 
quasi-equilibrium formulation. 

Before concluding this Section, it is important to note that 
the choice of a critical value for the axis ratio (a2/ai)crit also 
allows us to determine the ratio between the NS self-gravity 
and the tidal forces. Using a well-known Newtonian argu- 
ment, when the binary is at the separation rtide, the ratio be- 
tween the tidal and the NS self-gravitational force for a fluid 
element on the stellar surface when the binary is at the sepa- 
ration rtide is 

Mbh / ai " ' 

MnS V''tide 



= 7^. 



(36) 



so that when TZ — 1 the tidal and gravitational forces are 
equal. Using the affine model and considering the tidal ra- 
dius as the one at which a^jai — (a2/oi)crit = 0.44, we 
can compute the values of the two forces at tidal disruption. 
Doing so for the binaries considered in Table [T] we find that 



n ~ 0.59-0.70 for the T = 2 cases and 7^ ~ 0.46 - 0.47 for 
the r — 2.1h cases. Our tuning thus reveals that the tidal dis- 
ruption begins earlier than one would naively think and when 
the tidal force is only ~ 1/2 — 1/3 the self-gravitational one. 
The tidal force when the binary is at the separation r^ide may 
also be compared to the self-gravitational force of the star at 
isolation, i.e. when it is a sphere of radius i?NS- This amounts 
to calculating the ratio TV when a\ — > i?NS, '-e. 



Mbh /^-Rns 

. ''tide 



Mns 



7^^ 



(37) 



and enables u s to compare our results with those of Taniguchi 
et al.| ( |2008| l, where it was found that 7^' ~ 0.07. More 
specifically, for the binaries considered in Table [T] we find 
TV ~ 0.08 — 0.11, which is in good agreement with the afore- 
mentioned result. Stated differently, this reveals that the tidal 
disruption begins when the tidal force is roughly only ^ 1/10 
of the NS self-gravity at infinite separation. 

4. RESULTS 

Having tuned and validated the model, we will next con- 
sider its predictions for the baryon mass of the torus as a 
function of the mass ratio q, the stellar compactness C, and 
the BH spin a. Because this space of parameters is three- 
dimensional, it is more convenient to consider constant-spin 
slices and hence we will first comment on a fiducial case of a 
spinning BH with a = 0.4 and then discuss how these results 
change across the possible values of the spin. 

Most of our results are summarized in Fig. l3] which shows 
the baryonic mass of the torus in units of the stellar mass, 
Mft tor/-^^6.NS, as a function of the stellar compactness C and 
of the binary mass ratio q, with the data referring to a binary 
in which the BH has a dimensionless spin parameter a = 0.4. 
Quite clearly, the final mass in the torus varies considerably 
across the possible space of parameters and is systematically 
larger the smaller C. This is rather obvious: the smaller the 
compactness, the more "Newtonian" the star will be and thus 
with a smaller effective gravity at the surface. In turn, this 
means that, all else equal, it will be easier to disrupt it even at 



l^b.to/l^b.Ns 'or a=0-4 



l^b,tor/l^b,NS 

0.5 




Fig. 3. — Baryonic torus mass in units of the stellar mass Mtior/Mj, ns 
shown as a function of the compactness C and of the mass ratio g, for a I3H 
with spin parameter a = 0.4. 



large distances from the BH (i.e. 7-tide is comparatively large) 
and hence to produce a larger torus. 

At the same time. Fig. [3] shows that the mass in the torus 
will be larger when the BH and the NS have comparable 
masses. Also this result is quite obvious: the smaller the mass 
ratio, the more unlikely it will be for the star to be tidally dis- 
rupted and to be accreted "whole" by the BH. Putting things 
together, a BH-NS system with large mass ratio and small 
compactness maximizes the yields in terms of torus mass. 
For the same reasons, a binary with with small mass ratio and 
large compactness will yield the smallest tori. To fix the ideas: 
for a BH-NS system with a ~ 0.4, the toy model predicts that 
Mfc,tor/^4,NS ^ 0.5 when C == 0.10 and q = 0.33, while 
essentially no tori are produced for q < 0.14 and C > 0.14 
(cf. third panel of Fig. |4] where this data is also shown with 
contour plots). Overalfour toy model suggests that at least 
statistically a BH with spin larger than ~ 0.4 is necessary to 
produce any astrophysically relevant torus. 

The generic predictions of the toy model for a = 0.4 remain 
unchanged when considering also other BH spins, extending 
smoothly from smaller to larger spins. This is summarized 
in Fig. H] whose different panels from top to bottom refer to 
a = OlJ, 0.2, 0.4, 0.6, 0.8, 0.85, respectively. The baryonic 
mass of the torus is still in units of the NS mass and is re- 
ported as a function of the compactness and of the mass ra- 
tio, but it is shown through contour plots to better quantify 
the results. The numerical values of some representative con- 
tour lines are shown and allow for a direct measurement (the 
contours are equally spaced in a linear scale), while the thick 
and black solid line shows the area below which no torus is 
created {i.e. the "no-torus" area with A/f,_toi7^fh,NS < 0.01). 
Finally, shown with a horizontal dot-dashed line is the most 
likely mass ratio for a canonical 1.4 Af© NS. 

Moving from the top to the bottom of Fig. |4]it is easy to 
recognize that the maximum mass attained at the smallest 
compactness increases significantly with the BH spin, rang- 
ing from Mb,tor ^ 0.18Mb.Ns ^ 0.25 A/q for a = 0.0, to 
A^b.tor > 0.95 Afh^NS - 1.33 A/q for a = 0.85. At the same 
time, the "no-torus" area decreases and virtually disappears 
for a > 0.6. Stated differently, for sufficiently large BH spins 
a torus is always produced and with non-negligible mass. 
As an example, taking as fiducial compactness the canonical 
value of C ~ 0.145, the torus mass at the fiducial mass ratio 
goes from Mb tor ^ 0.06 AT;, ns ^ 0.08 A/© for a = 0.4, to 
ATb^tor ^ 0.24Mf,,NS ^ O.MMq for a = 0.85. 

A complementary view, because it refers to a different slic- 
ing of the space of parameters, is illustrated in Fig. |5] which 
is the same as in Fig. [4] but it shows the baryonic torus mass 
as a function of the BH spin and of the mass ratio for a 
fixed compactness C = 0.145 (left panel), or as a function 
of the BH spin and of the compactness for a fixed mass ra- 
tio q = 0.14 (right panel). Both panels of the figure are 
rather self-explanatory and underline what has already been 
discussed above: large tori masses are possible for BHs which 
are spinning sufficiently rapidly or for neutron stars which are 
not very compact (favoring stiff equations of state). 

In summary: considering an astrophysically realistic mass 
ratio q ~ 0.14 and a conservative value of the stellar com- 
pactness C ~ 0.145 (we recall that eve n larger values were 
recently suggested in Ozel et al. (2010 1), the predictions of 



the toy model are that the torus mass should be 



for BH spins < a < 0.85. Such masses are comparable 
but also smaller than the ones predicted by the analysis of 
unequal-mass NS-NS mergers carried out by Rezzolla et al. 
(2010) . 

5. AN INTUITIVE INTERPRETATION 

In the previous Section we have shown that the complex dy- 
namics of the tidal disruption and subsequent accretion onto 
the BH is well-captured by the simple assumptions needed to 
build our toy model. In what follows we will show that an 
even simpler framework can be built to explain at least quali- 
tatively the results of the toy model. 

We have already noted that binaries with less compact NSs 
produce bigger tori as these are more "Newtonian" and hence 
can sustain smaller tidal forces before being disrupted. Stated 
differently, less compact stars have larger tidal radii rtide- In 
the usual arguments this quantity is generally compared to the 
ISCO, and the standard line of arguments says that a BH-NS 
binary will produce a torus if rtide ^ Hsco- This reasoning, 
however, is inadequate for systems which yield low mass tori. 
An obvious failure of the argument is offered by a NS that is 
disrupted by its BH companion exactly at j'tide = Hsco- In 
this case half of the star would still be outside the ISCO and 
thus potentially capable of producing a torus. 

The necessary, but not sufficient, condition for a BH-NS 
binary to yield a torus is thus better expressed as 



rtide + Ql(nide) - nSCO 
2i?NS 



1 



ftide - nSCO 



2R 



>0. 



(39) 



•NS 



where the second expression is obtained after recognizing 
that at tidal disruption ai(7-tide) — 2i?NS- Expression ( (39] l 
has three important properties: it is dimensionless, it com- 
bines the three fundamental lengthscales of our system, and 
it essentially measures how many NS diameters fit between 
rude + oi (rtide) and rjsco (cf- Fig. [T]i. In other words, ([39]l 
quantifies "how much useful space" there is for the NS to form 
a torus after it is tidally disrupted. 

At this point it is natural to associate this quantity directly 
to the mass of the torus as expressed in units of the NS mass 



M, 



b,tor 



Ml 



b,NS 



1 



rtide - nsco 



2R 



'NS 



(40) 



A/f,,tor < 0.24Aff,_NS - 0.34 Afe 



: 



(38) 



where the exact proportionality will depend (albeit weakly) 
on q and C. Not surprisingly, ( |40l ) reproduces, at least qual- 
itatively, all of the phenomenology discussed before and pre- 
dicted by our toy model. As an example, for a fixed BH spin, 
and hence fixed risco, the torus mass will increase for less 
compact NSs since for these rtide will grow and more rapidly 
than i?NS- Similarly, for a fixed compactness, and hence for a 
fixed i?NS, the torus mass will grow with the BH spin as does 
the difference rtide ~ ^isco (Hsco decreases more rapidly than 

^■tide)- 

6. CONCLUDING REMARKS 

The production of a massive torus orbiting stably around a 
rotating BH is a necessary ingredient in all models that ex- 
plain SGRBs in terms of the coalescence of binary systems 
composed of a BH and a NS or of two NSs. The accurate cal- 
culation of this mass inevitably requires the use of numerical- 
relativity simulations, which however are still very expensive 
and have so far been applied only to a tiny patch of the pos- 
sible space of parameters. In the case of BH-NS binaries es- 
pecially, the space of parameters is particularly extended as it 
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Fig. 4. — Baryonic torus ma.ss in units of the stellar mass M;, tor/A:f6,NS shown as a function of the compactness C and of the mass ratio q. From top to bottom 
the different panels refer to different values of the BH spin (a = 0.0, 0.2, 0.4, 0.6, 0.8, 0.85) and the numbers on the iso-mass contours indicate the constant 
spacing. Shown with a thick and black solid line is the area below which no torus is created (i.e. the "no-torus" region with Mi, toi/A:?6.NS < 0.01), while shown 
with a horizontal dot-dashed line is the most-likely mass ratio for a canonical 1.4 Mq NS. 
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Fig. 5. — The same as in Fig.plbut considering different slices of the space of parameters. In pailicular, the left panel shows the baryonic torus mass in units 
of the stellar mass A/f, ,„, /A/;, ns as a function of the BH spin and mass ratio for C = 0.145, while the right panel shows baryonic torus mass as a function of 
BH spin and compactness for q = 0.14. 



involves the mass ratio of the binary g, the stellar compactness 
C, and the BH spin a. As a result, we presently have only a 
very limited idea of what are the likely torus masses that this 
process will yield and hence have a rather limited ability to 
assess whether or not the merger of a BH-NS system under 
astrophysically realistic conditions will serve as a robust sce- 
nario for the powering of SGRBs. 

To compensate for this lack of knowledge, we have devel- 
oped a toy model that allows for the computation of the mass 
of the torus without having to perform a numerical-relativity 
simulation. In essence, we model the NS in the binary as a 
tri-axial ellipsoid which is tidally distorted as it orbits in the 
tidal field of a rotating BH and as described by the relativis- 
tic affine model. When the star is disrupted, we decompose it 
into a system of non-interacting "fluid particles" which move 
on the corresponding geodesies. We therefore compute the 
mass of the torus as the integral of the masses of the particles 
which do not fall into the BH. The only free parameter in our 
model is the radius at which the tidal disruption takes place 
and which we tune in terms of the ratio of the semi-major 
axes on the equatorial plane and with the aid of numerical- 
relativity simulations. The tuning requires care, but allows us 
to reproduce with precision the majority of the data, some of 
which shows inconsistencies of their own. 

As it is natural for a semi-analytic approach, the model has 
a limited range of validity, which we have decided to set in 
the following ranges for the mass ratio, compactness and BH 
spin: 0.10 <q< 0.33, 0.1 < C < 0.16, 0.0 < a < 0.85, 
respectively. Overall, the toy model predicts that high BH 
spins, small mass ratios and small NS compactnesses all en- 
hance the mass of the remnant torus. As a result, tori with 
masses as large as Af(, tor — 1-33 Mq are predicted for a very 



extended star with a stellar compactness C ~ 0.1 inspiralling 
around a BH with dimensionless spin a — 0.85 and mass 
ratio q ~ 0.3. However, when considering a more astrophys- 
ically reasonable mass ratio g ~ 0.14 and a conservative but 
realistic value of the compactness C ~ 0.145, the predic- 
tions of the toy model set a considerably smaller upper limit 

ofAf6.tor<O.34Af0. 

All of the phenomenology discussed above has a rather 
intuitive interpretation and it is easy to show that the torus 
mass is directly related to how much of the star falls be- 
tween the tidal radius augmented of the NS semi-major 
axis and the ISCO. Hence, collecting the three fundamental 
lengthscales appearing in the system, the simple expression 
Mfcaor/A/b^NS oc [1 + (rtide - nsco)/2i?Ns] is able to cap- 
ture qualitatively the predictions that our toy model can make 
quantitatively. 

Toy models are by construction approximate representa- 
tions of much more complex phenomena and their predictions 
are therefore intrinsically accompanied by errors. Bearing this 
in mind, the toy model presented here can be further improved 
as new and more accurate numerical-relativity simulations are 
performed and as their level of realism increases with the in- 
clusion of microphysical EOSs, magnetic fields and radiative 
transfer This will represent the focus of our future work. 
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APPENDIX 
THE PARALLEL PROPAGATED TETRAD 



In this Appendix, we gather together the equations derived by Marck ( 1983 ) to define a tetrad which is parallel-transported as 
it moves along a timelike geodesic of a Kerr BH spacetime in Boyer-Lindquist coordinates 



ds^ 



1 



2MBHr 



At' 



^^asin^ dAtAd) + — dr^ 

S A 



I]d02 + A sin2 

2—i 



(Al) 
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where 



S = r2 



~2 2 

a cos ( 



A = r^ + a2 _ 2MBHr , 
A={r'^ + a^) - Aa^sin' 



(A2) 
(A3) 
(A4) 



In our model, the test particle is identified with the centre of mass of a NS orbiting its rotating BH companion. Marck expresses 
his results in the canonical symmetric orthonormal tetrad introduced by Carter (1968) 



u;W = 



u;W = 



,(2) 



A 

'A 

VSd6i, 



{At - a sin^ ed<j)) , 
dr , 



^(3)^smu^^^^ /„2 



[adt - {r^ + a^)d(j)] 



which has the convenience of casting the Kerr metric in the form 



ds^ = 



V{,.)i^ 



,a;('')a;M 



(A5) 

(A6) 
(A7) 
(A8) 



(A9) 



where il(tJ.){iy) — diag(— 1, 1, 1, 1) is the me tric tensor of Minkowksi spacetime. Before expressing the components of the basis 
vectors of the tetrad found in Marck ( 1983 i, we define the quantities a and /3 



K — a? cos^ I 



K 



K 



K — a? cos^ 



where K is Carter's constant, and the two vectors 



z (0) 



^M^^rr, 



, (1) ar[E{r^ + o?) — aL 



z (2) 



z (3) 



Vi^SA 
(3 a cos 9 [aE sin 9 — L^ sin^^ 



(AlO) 
(All) 

(A12) 
(A13) 
(A14) 



I3\l —acost 
K 



and 



_ (0) ar[E{r^ + d^) - dL 



VSA 



cr (1) 



Z (2) 



Z (3) 



dEsint 



(A15) 

(A16) 
(A17) 



L, sin ( 



where the dots indicate derivatives with respect to the proper time r, and where E and Lz are, respectively, the energy and the 
angular momentum about the axis of symmetry of the BH per unit mass of the star. We are now ready to express - in Carter's 
symmetric tetrad - the components of the vectors forming an orthonormal tetrad parallel-transported along an arbitrary timelike 
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geodesic in a Kerr spacetime. These are 



(0) _ E{r'^ + a^) ~ aL^ 



AVE 



(1) 



(2) 



(3) 




aE sin 6 — L^ sin 



(A18) 

(A19) 
(A20) 
(A21) 



ei = cos ^ei — sin ^62 , 
62 = sin 'I'ei + sin ^§2 , 



and 



(0) 



i^A 



I cos Or 



(1) a cos 9[E ir"^ + 0?) ~ aLz 



(2) 

63 = - 



(3) 



r{aEs\Ti9 — L^ sin^^ 1 



(A22) 
(A23) 



(A24) 
(A25) 
(A26) 
(A27) 



The rotation by an angle ^I* in (A22| and ( A23 1 ensures that gq) and 6(2) are indeed parallel-transported along any Kerr timelike 

itie 



geodesic. Finally, the evolution ot the angle m is governed by ([T9]l for circular equatorial orbits 
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